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The equilibrium magnetization of concentrated ferrofluids described by a system of polydisperse 
dipolar hard spheres is calculated as a function of the internal magnetic field using the Born-Mayer 
or cluster expansion technique. This paper extends the results of Phys. Rev. E 62, 6875 (2000) 
obtained for monodisperse ferrofluids. The magnetization is given as a power series expansion in 
two parameters related to the volume fraction and the coupling strength of the dipolar interaction, 
respectively. 
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I. INTRODUCTION 

Ferrofluids |^ are suspensions of ferromagnetic particles of about 10 nm diameter in a carrier fluid. The particles 
are stabilized against aggregation by coating with polymers or by electrostatic repulsion of charges brought on their 
surface. As long as the concentration of the particles is low, the equilibrium magnetization of a ferrofluid is that of 
an ideal paramagnetic gas. In highly concentrated ferrofluids on the other hand, the magnetization is influenced by 
effects of particle-particle interactions. 

We studied these effects for ferrofluids that are described by a system of identical dipolar hard spheres in , from 
now on referred to as paper I. In that paper we used the technique of the Born-Mayer or cluster expansion technique 
to evaluate the equilibrium magnetization as a series expansion in terms of the volume fraction (j) = NttD^ /6V, and 
a dipolar coupling parameter e = rri^ /Air^QkT, with N/V being the particle density, and D and m being the common 
hard sphere diameter and magnetic moment of the particles, respectively. 

However, real ferrofluids are polydisperse, i. e. the particles vary in size and magnetic moment. This property has a 
strong influence on the equilibrium magnetization, for concentrated as well as for dilute fluids. The goal of this paper 
is to generalize the findings of paper I to include the effects of polydispersity. 

The linear response problem of determining the static initial susceptibility x of a mixture of dipolar hard spheres 
was investigated already for the equivalent electric case in the framework of integral theories: The mean spherical 
model [3 w as extended to binary or multicomponent mixtures 0, | M 1 ^, IM 1^ - The reference hypernetted chain 
method |loj| was also applied to bi disp erse systems 0, . Recentlv |l3j| , the mean spherical model was used within 
the algebraic perturbation theory |l4| , however without leading to new results for the initial susceptibility. The mean 
spherical model was also extended to polydisperse ferrofluids in arbitrary high fields 0, |3 . Another theory dealing 
with arbitrary fields is the high temperature approximation |l7j| . A variant of this theory was proposed in |l8j | and 
extended in [l9l| . 

Our calculation follows closely that of paper I. Therein the application of the cluster expansion technique to a 
monodisperse system of dipolar hard spheres resulted in an expression for the magnetization M that can be put into 
the form 

oc 00 
m=0 r!=0 

where Mgat is the saturation magnetization of the fluid. The functions Lm,n{a) were given explicitly in terms of 
analytic expressions in the dimensionless magnetic field a. We calculated ^2,2(0:) and some of the Li_„(a). Lower 
orders vanish, except for the Langevin function Loyola). 

In the polydisperse case discussed here the parameters </>, e and a are replaced by more generally defined quantities 
(j), e and a (cf. Sec. Illlll . The calculated Lm.n transform into one-, two- or threefold sums over all particles, where 
the individual addends are analytical functions of the magnetic moments and diameters of the involved particles, and 
the reduced magnetic field a. 

The paper is organized as follows. In Sec. ^ we explain the principles of the cluster expansion technique. The 
main part of the paper is Sec. IIIII where we generalize the results of paper I for the equilibrium magnetization in the 
monodisperse case to polydisperse ferrofluids. The flndings are discussed in Sec. IIVI using example distributions. In 
Sec. the results are compared to experimental data. We conclude in Sec. IVII 
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II. CLUSTER EXPANSION: APPLICATION TO THE SYSTEM OF DIPOLAR HARD SPHERES 

Here we recapitulate briefly the principle of the Born-Mayer or cluster expansion technique: Consider a system of 
particles i = 1, ...TV interacting with an external potential Vi and with each other via a potential Vij. To calculate 
thermodynamic properties of the system one has to find the canonical partition function 

Z = [ e-^'='"=-^.<:' '''^ dT . (2.1) 



Here Vi — Vi/kT, Vij — Vij/kT, and dT means integration over the configuration space. The kinetic energy of the 
particles, if important, can be thought to be included in the terms Vi. One now writes 



where 



= e-^- - 1 . (2.3) 



If the typical interaction energy is small compared to kT, the fij can be considered as small parameters for the 
expansion of the integrand in Eq. H2.2|l . The leading terms factorize into low dimensional integrals that can be 
calculated at least numerically. 

In the system of dipolar hard spheres (monodisperse or polydisperse) the interaction potential Vij consists of a 
dipole-dipole (DD) interaction and a hard core (HC) repulsion part, Vij = Vj^^ + V^f^, where the first part is given 
by 

T.DD _ 3(m, ■r,j)(mj--f^j-)-m, -mj 



for two particles with magnetic moments m.i and m.j at a distance Vij = Xi — xj, with = |ry |, and fy = r^/ 



For particles with diameters Di and Dj one has V^'-^ {rij) = cxd, if < Dij = {Di + Dj)/2, and V^?'^(fy") = 



otherwise. 

Taking the thermodynamic limit in a system of dipolar particles requires some care because of the long range 
character of the forces .21] . We circumvented this problem by decomposing the dipolar potentials into a short range 
and a long range part, and replacing the latter by an effective mean field. Within this approach a particle experiences 
the local magnetic field 

M 

H/oca/ -\- H-dipole.near H ^~ ^" ^^dipole.near • (^■^) 

It consists of the dipolar near field 'H.dipoie,near that is produced by the other particles within a sphere of radius Rs 
and of an effective "external" field 

M 

H, =H+— , (2.6) 

seen by the particle in question at the center of the sphere. Here H is the macroscopic internal magnetic field and M 
the sought after equilibrium magnetization. Thus, when evaluating the partition function one has to take 

V = -m, • H, (2.7) 

as the external potential. The radius Rg of the sphere has to be taken to be sufficiently large to allow the far-field 
dipolar contributions to be replaced by those of a continuum - cf. paper I for details. Neither the kinetic energy of 
the magnetic particles, nor the carrier fluid has to be taken into account in the partition function, since these terms 
do not contribute to the equilibrium magnetization. The configuration space is thus given by the positions Xi of all 
particles and the orientations fii of their magnetic moments: dT = d'^'x.id^Vli. 

In paper I we used in addition also an expansion in the dipolar interaction: The /-terms were expanded as 

i^.=e-"re-"S°-l = i^^)+i^/)+/f) + ..., (2.8) 
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with 

= - 1 (2.9) 

= LJ!i2_Le--r ... ,n>l . (2.10) 
■' nl 

The two expansions concerning the /-terms and v^j^ together translate in the monodisperse case into a double power 
expansion of Z in the the volume fraction (p of the particles and the dipolar coupHng constant e. We calculated the 
terms in O(0e") and in 0{(j)'^e'^) of Z and from that the equilibrium magnetization in the same order. 

III. CALCULATING THE EQUILIBRIUM MAGNETIZATION 

A. Notation 

Consider a system of A'' spherical hard particles with diameters Di, Dj^ carrying permanent magnetic moments 
mi, mjv contained in a volume V and subjected to a magnetic field H. Let D and m be some "typical" values for 
diameters and magnetic moments that arc discussed further below. We then define the parameter (j) related to the 
volume fraction of the hard spheres and the dipolar coupling parameter e as 

- NnD^ 

^=-^ ' ^= =3 • (3-1 

The equilibrium magnetization is calculated as a power expansion in these two parameters. The dimensionless 
magnetic fields are defined as 



niH mHs 



kT ' ° kT 

Diameters and magnetic moments will be expressed in units of the typical values via 



(3.2) 



We will also use the minimal possible distance between two hard spheres i and j given by 



X,^liX + A,) = ^ . (3.4) 

Furthermore we introduce the reduced magnetic fields for each particle i by 

rriiH _ rriiH, _ 

oii = = Hia , Usi = = iJ-iUs . (3.5) 

Our cluster expansion docs not depend on how the characteristic values of D and m are defined in detail. For 
example they could be taken as as some weighted mean of the Di and r«,;, respectively, or their most probable 
values. To preserve this freedom of choice in our expansion oflfers some advantages for the comparison with magneto- 
granulometric analyses where the distribution of the diameters and magnetic moments is not known a priori but on 
the contrary the goal of the calculations. 

Note, however, that coincides with the actual volume fraction </> of the hard spheres only if one defines D via the 
mean volume of the particles 

^' = jjl2^i = [ ^'^(^) dD^{D^)^ . (3.6) 

i ^ 

Here P{D) is the normalized distribution function of the hard sphere diameters. 

Similarly m is related to the saturation magnetization Mgat of the ferrofluid via Mgat = Nm/fioV only if m is 
defined by 



m 



= ^ E = / m{D)P{D) dD . (3.7) 



February 6, 2008 



4 



The second equality of Eq. H3.7() holds when the magnetization of each particle is given by a function of its volume. 
We assume that this is the case and thus describe in this paper polydispersity effects of the ferrofluid by a distribution 
function P{D) depending only on the hard sphere diameter D. The generalization to a distribution function P{D, m) 
of independently varying diameters and moments is straightforward. Averages weighted with the distribution function 
P{D) of the diameters will mostly appear in the reduced version as integrals over the reduced diameter A = D/D 
with the appropriate weight function -P(A). 

The thermodynamic mean with respect to the noninteracting system will be denoted by (...)q and the corresponding 
canonical partition function by Zq. With this notation integrals over the /-terms appearing in Z ]2.2\i can be written 
in the form 

[l[e-^-f...dr^Zo{f...)o . (3.8) 

But in contrast to paper 1 we derive here an approximation directly for the free energy F = — fcTln^. If the 
particle-particle interaction would depend only on interparticle distance then F would be given by 

^ = ^^o-A:tQ^' + , (3.9) 

2 

including orders up to O(0 ), or, more generally speaking, up to terms of second order in the number density. The 
primed sums are taken over all particle pairs resp. all triples i,j,k. While Eq. H3.9|l does not hold for a system of 
dipolar particles in a magnetic field in arbitrary order of e it still is correct in the orders we want to calculate. 

The polydisperse generalization affects the calculation of the integrals in Eq. H3.9|l in two ways: (i) The fact that 
the individual dimensionless magnetic fields are different leads to more complicated expressions for some resulting 
functions compared to the monodisperse case - see the definitions of and below, (ii) The dispersion in the 
hard sphere diameters requires more difficult geometrical considerations concerning the w^'^-terms, especially in the 
three particle integral. 

B. The leading term: polydisperse Weiss model 

The leading term in Z is the partition function of the (formally) noninteracting paramagnetic gas in the magnetic 
field Hs 




(3.10) 



Zk = I e-'"'dxkdnk = 47rV^^^^^^ . (3.11) 

ask 



The equilibrium magnetization M{as) obtained from 



, 1 dF 1 d(kT\nZ) m d\nZ 



reads in leading order 



Here 



Nrfi 

M{as) ^ —^Cpoiyias) . (3.13) 
Cpoiyias) (^»"-) = / '"(^^'^ [M(A)a,] P(A)dA . (3.14) 

i 

is given by the sum of the Langevin paramagnetic contributions coming from each (reduced) magnetic moment 
fii — m.i/rn with C being the Langevin function. The second equality in Eq. H3.14|) is the continuous analog of the 
sum with /i(A) = m(A)/m and A = D/D. If one defines m via Nm — '^'^ that Nm/ iiqV = Mgat then 

'Cpoly(as ^ oo) = 1. 



February 6, 2008 



5 



The result H3.13|l reduces to the weU known expression for the magnetization of a polydisperse ideal paramagnetic 
gas as a superposition of Langevin functions, if one replaces ctg by ct (see e. g. jlfl]). However, the dipolar far-field 
contributions enter via H2.6|l as a mean field into 

as ^a+mM/3kT . (3.15) 

Thus, the lowest order result (|3.13() for M{a) 

Nm „ ( rr>,M\ 

M~/:po,.(^+^j (3.16) 

contains already corrections from the particle-particle interaction in the mean field approximation and (|3.16|l is the 
polydisperse generalization of the Weiss model [20l |. 

By replacing M on the right hand side of Eq. H3.16|l by the expression for the ideal paramagnetic gas, one arrives 
at the equation proposed in |l8j . 



C. The magnetization in 0{(j)) 

To calculate the canonical partition function in linear order of 4> we follow the lines of paper I, Sec. IV. One needs 
to include only the linear /-terms in the expansion H3.9|) . Thus we write 

= Fo - fcri^' . (3.17) 

For the second term in (|3.17|) , the trivial integrations over the degrees of freedom of all particles except i and j are 
performed first. This gives 

^E' (/'^)o = ^E'/ n^""/^^ dr (3.18) 

k 

= / e~^^^"^ fij dxidxjdflidQj 



Now we expand fij . Let 



such that 



An,^J = I e-^^-^^f^f dxdxjdnm, , (3.19) 



F-F,-Y,-—Y,A^,, . (3.20) 

n— 

We need not calculate since this term does not contribute to the equilibrium magnetization. Furthermore 
Ai = 0, because a dipolar magnetic field vanishes when averaged over a spherical surface. This is explained in more 
detail in paper I. Using the definition H2.10|l of fl^\ and the dipolar potential H2.4|l . we can write 

V f 
n\ J 

xe^"i2 r^jdrijdujijdflidrtj . 

Here we have integrated over Xi and decomposed r^j into the distance and a spherical angle w^ . The angles 
{ipi,'di), {ipj^dj), and (95,1?) represent the spherical angles fii, flj, and ojy , respectively. The function 

F((^j, Lpj,'dj,ip, d) = 3{m, ■ r,j){mj ■ r,j) - rh,, • rhj (3.22) 
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comes from the dipolar interaction. The integration over the directions of Vij , mi, and nij can still be done analytically. 
But in contrast to the monodisperse calculation, the result is now a function of two parameters a^i and agj. We 
define 



n\{n — 1)7T ZiZj 

X / g*^^^ COS l^j-j-Ctsj COS 19 j 

xP"((^i, t?j, (pi, t?j, (p, -d) dVLid^ljduJi 



(3.23) 



Gn{(Xsi, ctsj) is symmetric in its two arguments and a polydisperse counterpart to the function G„(q;s) defined in 
paper I. It is G^(q!s, ctg) — G„(q!s). Some of the are given in the appendix. 

Inserting Eq. 13.2311 into Eq. H3.21|l . integrating over between the minimal distance Z?y and oo, and introducing 
(j) and e yields 



and together with (|3.20|l the free energy 



F 



(3.24) 



(3.25) 



hT — I 

TT E '^^"E (M.M,)"A3-3«G^(a,„ a,,) 



N 



Here the dots represent the contribution from Aq that was not calculated. It can easily be shown that F does not 
depend on a particular definition of to or D. 

Now, the equilibrium magnetization M{as) is given in 0{(l)) by 



The function G'^^y^ „ is the derivative of 



Nrn 
1^ 



Cpoiy{as) + ^ 0e"Gpoiy,„(as) 



n=2 



(3.26) 



Gpoiy,n(as) = (AijAij)"Af^- ^"G^(as,,Q;s 



[M(A0A*(A,)]"A3r3n 



X G^(a,„a,,)P(A,)P(A,)dA,dA, 



(3.27) 



which is is a generalization of G„ T| and reduces to the latter in the monodisperse case Di = D = D and nii ~ m ~ m. 

In a last step, we convert the expression for M as a function of into a function of a using the definition of Hg 
in Eq. (|2.6(l . By expanding and iterating in a way that is analogous to the procedure in paper I, Sec. IV C we obtain 
the final result for M up to order (j) 



M{a) 



Nm r 



'Cpoiy(a) + 80e£poiy(a)£poiy(a) 

oo 

+ 5^^e"G;„iy,„(c 



(3.28) 



The two leading terms can be seen as the polydisperse extension of the high temperature approximation derived in 
[TtI for monodisperse systems. 



D. The contribution in 0{(f> 



For the monodisperse system, the magnetization contribution in 0(0 e^) was calculated in Sec. V and appendix B 
of paper I. The cluster integrals needed in that order are shown in Fig. 3 of paper I. Some of them vanish for the 
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same reason as the contribution Ai in 0{(j)€): They involve the averaging of a dipolar field over a spherical surface. 
Most of the remaining integrals cancel when the free energy F = — kT In Z is calculated. Up to O(e^) we can write 
the remaining term in Eq. 1)3. 9|l as 

If f f \ - / AO) AO) Ao)\ , o / f(2) .(0) .(o)\ 

\JijJikJjk/Q — \^Jij Jjk /f^' yij Jik Jjk 

+3(4/^i^^/;^')^ . (3.29) 

These three terms correspond to the graphs E, G and H in paper I. The first one is an 0{(j> e°)-term that does not 
contribute to the magnetization. 



1. Graph G 

Here we give an outline of the polydisperse generalization of the calculation pertaining to graph G (appendix B.7 
of paper I). The quantity Zq calculated in paper I, is given by the sum 

Zo-zJ^{f^ff^f^)^ (3.30) 

over distinct particles fc. For calculating i^fij'' ffk f^i^ O'^^ starts with the trivial integrations: the degrees of 

freedom of all particles except z, j, and fc, the position of the center of mass of the three remaining particles, and the 
orientation of rtifc, since particle k is not involved in dipolar interactions in this cluster. Then, for a fixed distance 
the integrations over Wy , Vti, and 17^, defined as in Sec. IIII Cl are carried out. This introduces the function Gf into 
the result: 

36 , , n— 2 



/ .(2) „(0) .(0)\ 30 2_2 p, , 

X / r^er'-^f^l'S dr^.dv,, . (3.31) 



The integral over r^fc can be described by the following geometrical considerations: The volume of possible positions 
of particle k has to be found, such that this particle overlaps with both, particle i (i. e. rik < Dik) and j (rjk < Djk). 

Otherwise the integral would vanish because of the factor f'l^ ■ This is only possible, if Vij < Dik + Djk- 
In a final step the integration over between Dij and Dik + Djk is carried out. The final result is 

x/«(A,,,A,fc,A,fe) . (3.32) 

The function f^{Aij, Aik, Ajk) is given in the appendix. The contribution to the free energy is according to Eq. H3.9|) 
and Eq. (|3.29|l given by 



/o 

= -NkTji' ^ ' 7''^ Gpoiy.2(5s) . (3.33) 



^_2_2H-61n2; 



The function 



Gpo.„2(S.) = (3.34) 

was defined in such a way, that it reduces to G2(as) in the monodisperse case. Finally introducing the diameter 
distribution function P(A) requires the replacement 

^ j F(A,) dA,P(A,) dA,P(A,) dAk (3.35) 

in eqs. H3.33I3.34|I . 
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2. Graph H 

The integrations to calculate \ fiP fjV fj^/ are performed as follows (compare to appendix B.8 in paper I): After 

doing the trivial integrations (concerning the possible configurations of the particles except i, j, and fc, and the center 
of mass of the cluster), the possible orientations of the cluster are integrated out. Then, the integrations over the 
orientations of mi, nij, and are performed. One arrives at 



/ Ai) AD A0)\ 2 -2 47r2 „ 



^(3cos^.,.-l)e-^re-"S^/(^) 

' ij ' ik 

xdr^jdr.ksm^jkddjk ■ (3.36) 

This is up to a factor N'^Zq/2 a polydisperse generalization of Eq. (B21) of paper I. The remaining integrations 
concern the distances r,j and r^^., and the angle "djk between r^j and r^^. The function is defined by 



K {asi, asj,ask) = 



8 ZiZjZk 

1 /.I 

X (3 + uf) UjUk duidujduk 



-1 j-i j-i 



[compare to Eq. (B22) of paper I]. An expression for the monodisperse counterpart K(as) was given in appendix A 
of paper I. Here we express via Langevin functions: 

K^{asi,asj,ask) = 

3C{asj)C{ask) [3 + /:'(as,:) + ^as^ f] . (3.38) 



Note that > Dij and r^fc > Dik is required, otherwise the integral (|3.36() vanishes. The requirement that particles 
j and k have to overlap imposes the additional restrictions \rij — rik\ < Djk and 

^jk < ^fk" = arccos '^J'^ ~ ^'"^ . (3.39) 

After performing the remaining integrations within these limits the result is 

48 



/ADfWA^)\ _ „2„ „ T-2-2 

\hj hk Jjk - e 



5A^2 

xiC^(a.^,a,j,a.fc)/^(A»j,A,fc,Aj7c) • (3.40) 
The function /^'(A^j, A^^^, Ajk) is given in the appendix. The contribution to the free energy is 

(fS^fikfl^)^ = iVfcr0Vifpoiy(a.) , (3.41) 

with 

Kpo\y{as) = -y^X! (3.42) 
xii'^(as,, a^j-, asfc)/^(Ay, Ajfc, A^-^) . 



3. The magnetization contribution 



Inserting Eq. I|3.33|l and Eq. (|3.41|) in Eq. H3.9|l and calculating the equilibrium magnetization results in an additional 
O(0^e2)-term in Eq. IpHH]) : 



Nm—2_ 
m 6 



1 + 61112 ~, 



(3.43) 



February 6, 2008 



9 



The expansion and iteration procedure to switch from as to a is identical to the monodisperse case in paper I. The 
fuU expression for the magnetization containing all calculated terms reads 



Nm 
M = — 



'Cpoiy(a) + 8(/)e£poiy(a)>Cpoiy(a) 

oo 
n=2 

+Q0e^C^oiy{ol)C'^,,y{af (3.44) 
+32(/) e^£poiy(a)^£poiy(") 

Gpoly,2(a) - T'^'^ K'^o\y{0l) 



IV. SELECTED RESULTS FOR THE MAGNETIZATION 

The figures ^ and [21 showing separate contributions to M{a) in polydisperse systems were obtained by using 
lognormal distributions 

P(D) = e--"/^e~ i"'(^/A,)/(2-^) (4.1) 

for the particle diameters as a representative and often used example for size distibutions of model polydisperse 
systems. Using similar distributions like, e. g., Gamma distributions, however, would not qualitatively modify the 
results. In Sec. we use also an experimentally determined size distribution. Here, the quantity D was taken to be 
defined via the mean volume (|3.6|l . i.e., D = (L)'^)^ = D^e^^'^ so that (j) is the volume fraction 0. The magnetic 
moments of the particles were taken to scale with their volumes, m ~ D^, allowing to set /i(A) — such that 
TO — {m)p. 

Figures n and 12 show the contributions 

Li,i(a) = 8/:poiy(a)/:;„iy(s) , (4.2) 

^2,2(0) = 64/:poiy(a)£poiy(a)2 

+32/:poiy(a)2/:;,'oiy(s) (4.3) 
l + 61n2 ~, , 

H t^poly,2(a) - ■^polyi") . 

ii,2(a) - G'p^iy^^ia) , (4.4) 
ii,3(a) = G;,,iy_3(S) (4.5) 

to M (a) (|3.44|) for different values of the width a of the distribution l|4.1|l . The contributions of higher-order terms 
increase with growing a. This is so because they depend on higher moments of the distribution P(A) that grow with 
the width of the distribution, even if the third moment ^ A'^^^ is kept fixed. The shift of the maxima of the curves to 
the right has a similar reason: Bigger particles, that react to smaller fields get more and more important when the 
width of the distribution grows. 

The assumption /x(A) = A'^ that the magnetic moments in our ferrofluid model scale ~ with the total volume 
of its hard sphere constituents is somewhat too simple for particles in real ferrofluids for two reasons: First, in the 
common case of steric stabilization by polymers surfactants providing the repulsion the surfactant layer of about 1 
- 3 nm does not contribute to the magnetic moment. Second, an outer layer of the magnetic material might be 
magnetically dead so that it does not contribute to the magnetic moment either. For magnetite particles a dead layer 
depth of « 0.8 nm has been reported [S^. To account for the sum of these two effects we have introduced in our 
calculation for M an effective magnetic diameter -Dmag via the relation D — I?mag + 5.6 nm. It ascribes to every 
particle that has a magnetically effective core of diameter Dmag a hard sphere with a magnetically inert layer of depth 
2.8 nm. The magnetic moment of each particle was then taken to scale with is magnetically effective volume, i.e., 

/i= (^mag/^D)^. 

The full line in Fig. |3| shows the reduced equilibrium magnetization of an interacting polydisperse ferrofluid as 
function of a. Here each particle has a total magnetic inactive layer of 2.8 nm and the diameters -Dmag of the 
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magnetic cores are lognormally distributed with a = 0.25 and (-Dinag)p"^ ~ 10 The latter is taken as the 

reference diameter D in the calculation. The particle density is chosen in such a way that cj) = 0.05. Note that here 
= 'Amagj the volume fraction of the magnetically active material. The magnetic moment to for Dmag = 10 nm is 
chosen such that e = 2. The magnetization of the core is then about 550 kA/m at room temperature which is a little 
bit larger than that of magnetite. 

Expansion terms of order (/)e" are taken into account up to order n—5 for the magnetization curves in Fig.|31 Higher 
O(0e")-terms have only a small effect on the magnetization. Here e is by definition of m and D a typical interaction 
energy divided by kT for particles at a distance of Z3 = 10 nm. But with the two additional dead layers of total size 
5.6 nm in between the particles of our model ferrofluid the real typical dipolar energies at contact are smaller. 

The magnetization is compared to that of a polydisperse fluid without particle-particle interaction (long dashed 
line), and to a monodisperse ferrofluid, both with and without taking into account the particle-particle interaction 
(short dashed and dotted line, respectively). The monodisperse system consists of particles with Dmag ==10 nm with 
the same nonmagnetic layer thickness, bulk magnetization, and 0,nag as before. 

One sees that taking into account polydispersity or particle interaction alone strengthens the magnetization and 
especially the initial susceptibility. Both effects are comparable for the given parameters. Together, they result in an 
even higher equilibrium magnetization. 

Figure 0] shows magnetization curves for magnetite-based ferrofluids with a bulk magnetization of 480 kA/m for 
distributions of different widths, i^mag is taken to be lognormally distributed with (£'mag)p = (8 nm)'^ and cr — 0.2, 
0.3, and 0.4. The volume fraction of the magnetic material is 0mag = 0.1. The particles are again assumed to carry a 
nonmagnetic layer of 2.8 nm thickness. 

The increase in a causes an increase of the initial susceptibility already in the noninteracting case (long dashed 
lines). Including the 0((/)e)-terms (short dashed lines) has a positive effect on the magnetization. The relative increase 
is maximal for small a. The magnetization decreases again at higher a, if higher order terms are taken into account 

(solid lines). For the considered ferrofluids the 0((/) e^)-term that is negative for higher a (see Fig.^ is almost solely 
responsible for this decrease. The positive contributions from the higher O(0e")-terms (n > 2) are again negligible, 
except for cr = 0.4 and small a, where they cause a further increase of the initial susceptibility. For small a the 

O(0 e^)-term also has a positive effect, but this effect is too small to be visible. The plots show again, that the 
influence of higher order terms is larger for broad distributions. 

V. COMPARISON TO EXPERIMENTS 

We compared our theorerical predictions for the magnetization curves with experimental results of two different 
magnetite-based ferrofluids. S. Odenbach (ZARM, Bremen) provided data on the equilibrium magnetization of the 
ferrofluid EMC 905 produced by FerroTec. We fitted our theoretical result (|3.44|) taking into account the terms in 
0((/)e") up to n = 5 to the data assuming lognormally distributed Mnag and a nonmagnetic layer of depth 2.8 nm. 
The bulk magnetization of magnetite was taken to be 480 kA/m. The result is shown in Fig.O According to the fit, 
the saturation magnetization of the ferrofluid is Mgat — 37 A kA/m. The parameters defining the distribution turn 
out to be Dq = 8.3 nm and a = 0.28. There are small differences between the data and the fit curve that are in our 
opinion due to deviations of the real diameter distribution curve from the idealized lognormal form. 

J. Embs (Universitat des Saarlandes, Saarbriicken) measured the equilibrium magnetization curve of the ferrofluid 
APG 933 of FerroTec. In addition he determined the diameter distribution of its particles by transmission electron 
microscopy (TEM) which was then used in our theoretical analysis. Diameters found in TEM measurements are 
those of the magnetite particles. We assumed the magnetically effective diameters Z?niag to be 2 x 0.8 nm = 1.6 nm 
smaller and to be zero for particles smaller than 1.6 nm. The hard core diameters D were taken to be 2 x 2 nm — 4 
nm larger than the diameters obtained from the TEM measurements. As above, we took into account terms up to 
0((/)e^) and set the bulk magnetization of magnetite to 480 kA/m. Fig. |H|shows the TEM data and the experimental 
magnetization curve together with the results of our theory. Both agree very well. 

VI. CONCLUSION 

In this paper we used the technique of cluster expansion to derive an approximation to the equilibrium magnetization 
for the system of dipolar hard spheres in a magnetic field with diameter and/or magnetic moment dispersion as a 
model system for a polydisperse ferrofluid. The calculation results in an expression for the magnetization M in form 
of a twofold series expansion in the parameters 0, closely related to the volume fraction 0, and e, a coupling parameter 
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measuring the strength of the dipolar interaction: 

^^-^EEr^"w(«) . (6.1) 

4>, e, and the dimensionless magnetic field a are defined for some typical values D and m for the hard sphere 
diameters and magnetic moments respectively, m can be chosen in such a way that the prefactor Nm/fj,QV reduces 
to the saturation magnetization of the system. We gave expressions for (n < 5) and ^2,2- Lower orders vanish, 
except for Lq^q reducing to the Langevin function in the monodisperse case. The calculated im,n can be written as 
multiple sums over all particles whose addends are analytical expressions. 

The influence of particle-particle interaction grows with increasing width of the considered diameter distribution. 
Taking into account only the Li i-term results in an increase of the magnetization relative to the non interacting sys- 
tem, whereas the L2,2~term leads again to somewhat smaller values at higher a. Only at very small a its contribution 
is positive. The Li „-terms have little effect for realistic, magnetite-based ferrofluids, except for broad distributions, 
where they increase the initial magnetization. 
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APPENDIX A: THE FUNCTIONS G^, /^, /^ 

The functions G^{xi,X2) are symmetric in their arguments and have the form 

+ coth(a;i)G^W(l/2;i,l/a;2) (Al) 

+ coth(a;2)G^W (1/2^2, IM) 

+ coth(a;i)coth(a;2)G^(2)(l/a;i, 1/2^2) ■ 

The Gn^^'^ are polynomials and read for n < 5 

r<P'S>)( \ 8 , 4 2 4 2 12 2 2 
G2 ^ '{yi,y2) = - + -Vi + -2/2 + yyi2/2 , 

G^2^^\yuV2) = -72/1-^2/12/1 , (A2) 



G^'^ (2/1, 2/2) 



5" 5 
12 

"5" 



G3' '(2/1,^2) = -—2712/2 - —2/12/2 

-^2/12/2 - y 2/12/2 , 
4 24 

6^^^(2/1,2/2) - -^2/2 + ^2/^2/2 (A3) 

+ ^2/| + y2/i2/| , 

r^P(2), N 16 4 2 4 2 12 2 2 
G3 (2/1,2/2) = ^ - ^2/1 - ^2/2 - y 2/12/2 . 
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It is 



and 



PiO)r ^ 8 4 2,4 4 4 2 

4 (2/1,2/2) = + + + 



12 2 2 36 4 2 44 

+ y2/i2/2 + —2/12/2 + 352/2 
+y2/?2/2' + 122/^2/2^ , 
^(2/1,2/2) = -^2/1-352/1-352/12/2 

Gr'^(yi,y2) = ^2/12/2 + y2/?y2 + ^2/12/1 + 12j/?y^ 



.P(o). ^ 12 104 3 60 5 

^5 (2/1,2/2) = ^2/1^2 - ^2/1 2/2 - 2/1 2 



_ 104 3 _ 732 3 3 _ 240 5 3 

60 . 240 3 . 540 . . 
-yy2/i2/2 - ^2/12/2 - ^2/12/2 , 

^P(i)/ ^ 4 4 2 60 4 

^5^^2/1, 2/2) = -23^-2/2 + 3352/1 2/2 + y;^2/i2/2 

4 3 172 2 3 240 4 
+ 3552/2 + ^2/1 2/2 + ^2/1 2/2 

4 , 60 2 5 540 4 5 
+ y;^y2 + yy2'i^2 + ^2/1 2/2 , 

^■P(2)/ ^ 16 , 8 2 4 4 

8 2 ^2 2 2 4 2 

-j^ -j^ g ^ 2/2 rji-^y^^y^ -j^-j^2/i2/2 

44 60 2 4 540 4 4 
— —2/2 ~ Yj2/i2/2 jy2/i2/2 • 



2 



m,,A.,A,.)=ln(^^)-^(^^ 

8(Af, + A3,) 3(A^, + Ag,) 

3A3^ A^ 
13A^, + 14A,fcA,.fc + 13A|, 



2 



12(Aife + A,-fe) 



/ (A,j, Aife, A,fc) - -g - ^ ^3;^ 



A|, 2A3, 2A3, 



4A,,A,fc 9A3. 9A3, 
, (At,+A^ + Aj,)(Ag^+A^ + Ag,) 



8A?,A3, 



When all diameters are equal, these functions reduce to 



/G(A,A,A) = i+ln2 
6 

/^(A,A,A) = -A . 
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FIG. 1: The contributions Li,i and L2,2 as functions of a for lognormal distributions with different widths a. Here (A^)^ = 1 
and /j(A) = A^. 
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FIG. 3: Equilibrium magnetization of a polydisperse ferrofluid ((Omag)p = 10 nm, a = 0.25, n = (Dn^ag/O)-') and a 
comparable monodisperse ferrofluid, both with and without taking into account the particle interaction. See text for further 
details. 
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FIG. 4: Equilibrium magnetization as a function of a for lognormal distributions with (-Dmag)p ^ = 8 and different widths 
cr. Long dashed: without particle interaction, dashed: particle interaction only in O(0e), solid lines: all calculated terms. 



February 6, 2008 18 



40 



30 



E 

3 20 



10 








H (kA/m) 

FIG. 5: Equilibrium magnetization of the ferrofluid EMC 905. Dots: experiment, solid line: the theoretical magnetization 
curve assuming lognormally distributed Dmag, c.f. text. 
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FIG. 6: Top: equilibrium magnetization of the ferrofluid APG 933. Circles: experiment, solid line: the theoretical magnetization 
curve. Bottom: diameter distribution according to TEM measurements. 



